Realization of holonomic single-qubit operations 
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Universal single-qubit operations based on purely geometric phase factors in adiabatic processes 
are demonstrated by utilizing a four- level system in a trapped single ''"Ca"'' ion connected by three 
oscillating fields. Robustness against parameter variations is studied. The scheme demonstrated 
here can be employed as a building block for large-scale holonomic quantum computations, which 
may be useful for large qubit systems with statistical variations in system parameters. 
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I. INTRODUCTION 

Experimental studies of quantum inform.ation process- 
ing (QIP) has progressed much in recent years. In stud- 
ies of QIP using trapped ions, realization of small scale 
computation and entanglement generation with relatively 
high fidelity has been reported |1[, which reveals that 
there is no fundamental obstacle to scaling a large num- 
ber of ions. The upcoming challenges are large-scale op- 
erations and high-fidelity gate operations toward fault- 
tolerant quantum computation. 

As a way to realize high-fidelity gate operations, quan- 
tum gates and quantum computation using geometric 
phase factors have recently been studied. This orig- 
inates from Holonomic Quantum Computation (HQC) 
proposed by Zanardi et aL[2j In HQC, degenerate mul- 
tiple quantum states are utilized and unitary operations 
are performed by varying the system Hamiltonian along 
a closed path in the parameter space. The final state in 
the HQC is dependent only on the global property of the 
closed path; therefore, HQC is considered to be robust 
against certain types of errors. Even when diabatic evo- 
lutions of the system or nondegenerate quantum states 
are used, similar advantages can be expected as long as 
the unitary operations performed are determined by geo- 
metric phase factors [Geometric Quantum Computation 

(GQC)a]]- 

There have been a number of experiments in different 
systems related to GQC (JtMI- It is known that univer- 
sal quantum computation can be realized with a combi- 
nation of single-qubit and two-qubit operations [9]. As 
demonstrations of two qubit operations in trapped-ion 
systems, "Geometric-phase gate" by Leibried et al. [5| 
and the M0lmer-S0rensen gate [3, [Ifl, [III have been re- 
alized. A gate fidelity of 99.3% has been realized using 
such a scheme @. 

Single qubit operations in trapped-ion systems have 
been performed to date using a dynamical method with 
variable pulse lengths and phases. By replacing such a 
method with those that use geometric phases, gate op- 
erations that are robust against variation of parameters 
such as pulse intensity and lengths can be expected. 

In this work, we report the realization of purely geo- 
metric single qubit operations using a four-level system 



in a single ^°Ca+ ion. Three transitions in the four- level 
system are excited either by three optical fields or by two 
optical fields and one RF magnetic field. Rotation opera- 
tions by arbitrary angles along two different axes, X and 
Z in the Bloch sphere for the qubit, are demonstrated by 
utilizing stimulated Raman adiabatic passage (STIRAP) 
in this four-level system. Robustness against parameter 
variations is also demonstrated. 

This work is based on proposals using a tripod system 
comprising of one upper state and three lower states con- 
nected by three oscillating fields |12l4l5| . There are two 
dark states in the system, and those dark states are adi- 
abatically manipulated with the intensities and phases of 
the oscillating fields to perform single-qubit operations. 

There was also a proposal and demonstrations of 
single-qubit operations using geometric phase factors us- 
ing a two-level system and square pulses [y, Q- How- 
ever, to the best of our knowledge no adiabatic (or holo- 
nomic) demonstration of single-qubit operations using 
dark states has been reported so far. 



II. PRINCIPLES FOR HOLONOMIC 
SINGLE-QUBIT OPERATIONS 

Here the formalization of the scheme is summarized 
based on Kis and Renzoni [13] for later reference. A 
four-level system comprising of {|0) , |1) , |2) , \u)} is con- 
sidered [see Fig.[T][a)]. |0) and |1) span a qubit manifold, 
while \u) represents an upper state and |2) represents an 
auxiliary state used in STIRAP. Three near resonant os- 
cillating fields are applied to this system between \u) and 
|fc) (fc = 0,1,2), and the Hamiltonian in the interaction 
picture with the rotating wave approximation is given as 



H{t) - -hS \u) ("I + 2 E [^'^ W I") (^1 + H-C-] , (1) 

A:=0 

where S is the detuning of the lasers, which is assumed 
to be common to all the three transitions. The Rabi 
frequencies for the transitions between \u) and the qubit 
states are il.o{t) = U{t)cos6i (assumed to be real) and 
rii(i) = ri(i)e"^i sin 6*1, where a real quantity fl{t) is the 
common envelope function. 6i = tan~^ |f2i(t)/f2o(i)| and 
4>i = a,TgQi(t) are assumed to be independent of time. 



The Rabi frequency between \u) and |2) is described as 
n2{t) = \rt2{t)\e'*^<-'^ with 02 (i) = argf^aW- 

The fohowing states in the qubit manifold are use- 
ful for describing the Hamiltonian and explaining the 
gate procedure: |Ci) — cos0i|O) + e"*"^! sin6'i |1) and 
\Di) ^ - sin 6*1 |0)+e-*"*i cos 6*1 |1), and similarly, the fol- 
lowing states in the total lower-state manifold spanned by 
{|0) , |1) , |2>}: IC2) = cos02(i) \Ci) + e'^-^^W sin^2(i) |2) 
and \D2) ^ - sm02{t) |Ci) + e-*"^^(*) cos6'2(t) |2), where 
02 (i) — tan^^ |J72(i)/i^(i)|- Based on these, the Hamilto- 
nian can be simplified as follows: 



Hit) = ~h5\u) {u\ + -rftotalW i\u) {C2\ 



H.c. 



(2) 



,1/2 



where f^totaiW = [Q^t) + 1^2 {t)\^Y'^ . When 

S — 0, this Hamiltonian has four eigenvec- 
tors {\Di) ,\D2) ,\B^) ,\B_)} with eigenvalues of 
{0, 0,firitotai(t)/2,-firitotai(i)/2}, respectively, where 
|i?±) = (IC2) ± \u))/V2. 

We proof the ability to perform arbitrary single- 
qubit operations by first starting from an arbitrary ini- 
tial state in the qubit manifold, \tpo) = a |0) -I- /3|1), 
which can be rewritten using |Ci) and \Di) as |'0o) = 
(Ci|-0o) |Ci) + (Dilipo) \Di). This initial state wiU be 
transferred in the first STIRAP sequence that makes 
use of \D2} [the left half in Fig. ^h), with 6*2 (i) : 
tt/2 — ?► and (f)2 = argJ72 = 0] to the intermediate 
state IV-i) = - (CilV^o) |2) + {Di\ijjq) \Di). This inter- 
mediate state will then be transferred using the sec- 
ond STIRAP sequence [the right half in Fig. [Hb), with 
02{t) : — > 7r/2 and ^2 = a.vgVL2 — $] to the final 



state [■(/'2) 



(CilV'o) \Ci) + (T^ilV-o) \Di), which can 



be rewritten using the identity operator / and the Pauli 
operators a = {cTx, ^y,^z) on the computational manifold 
{[0) , |1)} as 1^2) = e**/2[cos($/2)/-|-m-asin($/2)] [V'o) 
with n = (sin 2^1 cos (/)i, — sin 26*1 sin 01, cos 26*1). This 
operation corresponds to a rotation by angle — $ around 
n, which can be taken arbitrarily by selecting the values 
for 01 and 4>i , hence arbitrary single-qubit operations can 
be performed. 

This scheme is based on adiabatic population trans- 
fer using I-D2), and diabatic transitions from this state 
to \B±) are among the possible causes of infidelity 
in gate operations. The probability for such dia- 
batic transitions is calculated to be in the order of 



202{t) /nl,,, 



[l6[ , which should be set to be much 



smaller than 1 to maintain high fidelity. 



III. EXPERIMENTAL SETUP AND 
PROCEDURES 

The experimental setup has been previously described 
[iTj and only a brief description is given here. A sin- 
gle ^°Ca+ is trapped in vacuum (6 x 10~^ Pa) using a 
linear Paul trap. The trap used here is a conventional 
linear trap with an operating frequency of 23 MHz and 



secular frequencies of {ujx,ujy,ujz)/2TT = (2.4,2.2,0.69) 
MHz. A bias magnetic field of 2.9 x 10~* T is applied 
to define a quantization axis, which results in a Zeeman 
splitting of ^4.9 MHz between I?5/2 sublevels. In ex- 
periments that employ three optical fields, a titanium 
sapphire laser at 729 nm stabilized to a high-finesse low- 
thermal-expansion cavity is used for the excitation of ions 
between 5*1/2 and D^/2- The amplitudes and frequencies 
of the three optical fields at 729 nm are changed by vary- 
ing the three RF fields that are combined and fed to an 
acousto-optic modulator. The RF fields are generated 
by three direct-digital synthesis (DDS) boards that are 
controlled by a field-programmable gate array. Polar- 
ization of the optical fields is adjusted so that there are 
polarizations both parallel and perpendicular to the bias 
magnetic field; therefore, the transition to be excited is 
selected by changing the frequencies of the DDS boards. 
In experiments that employ two optical fields and one 
RF magnetic field, an RF coil in the vicinity of the trap 
is used to generate the RF field. The details of the RF 
excitation procedure are similar to that described in Ref. 

m- 

We have chosen the Sii2~D^/2 electric-quadrupole 
transition of ''"Ca"'" for realizing the geometric phase gate. 
The encoding of the tripod system to the sublevels in 
5i/2 and -D5/2 is as follows: Si/2(jnj = —1/2) as the 
"upper" state \u) and three Zeeman sublevels in D^i2 
{rrij = —3/2,1/2,-5/2) as the lower states |0), |1) and 
1 2), respectively (see Fig. [2c)). The reason that a level 
in S'1/2 is adopted instead of one in P3/2 as \u) is to 
avoid the effects of spontaneous emissions when all the 
fields are resonant. However, unwanted couplings be- 
tween 5'i/2(toj = 1/2) and |0), |1) and |2) must still be 
considered. Since these couplings are off resonance, vary- 
ing 02{t) leads to time-dependent AC-Stark shifts that 
disturb the null eigenenergies of the qubit states. This 
effect was numerically evaluated with realistic parame- 
ters and the phase accumulated during the gate operation 
was confirmed to be smaller than 0.5 rad under typical 
experimental conditions in the present work. 

We essentially do not compensate such AC-Stark shifts 
for the present work. In the experiments in ^IVI and of 
rotations by angle n in i fVl we adjust one of the detunings 
(one for |u)-|0)) so that fringe shift resulting from AC- 
Stark shifts are apparently canceled. This cancelation 
depends on the details of the STIRAP pulses, including 
such conditions as peak Rabi frequencies, pulse shapes or 
total time (it may not depend on such conditions as input 
states or rotation angles since AC Stark shifts do not 
depend on these). We also evaluate gate performances 
without cancellation of this type by examining rotations 
by angle 7r/2 with relatively low Rabi frequencies in JfVl 
andflYnil 

The procedure for gate operations is as follows. Firstly, 
the ion is cooled to near the motional ground state with 
Doppler cooling by 397 and 866 nm and with sideband 
cooling of the axial motion by lasers at 729 and 854 nm. 
The average motional quantum number along axial direc- 



tion after sideband cooling is fiz = 0.06. Optical pumping 
using a 397 nm a~ transition is then performed to initial- 
ize the ion in Si/2{'mj — —1/2) (|m)). The ion is trans- 
ferred from \u) to the computational manifold spanned 
by {|0) , |1)} using square ir pulses, and then a gate STI- 
RAP pulse sequence is applied between |0) and |1). The 
pulse shape of the gate consists of partially overlapping 
sinusoidal curves and constant values, as in Fig. [Ijd). 
A step variation of 4'2it) by $, which brings a geometric 
phase factor, is given at the middle of the gate sequence 
where |^2(0I = 0- After the gate sequence, a square 
TT pulse is applied to map the state of |0) to |u). An 
additional pulse may also be applied, depending on the 
element of the density matrix to be observed. State dis- 
crimination is performed by detecting fluorescence from 
the ion with a photomultiplier tube during a period of 7 
ms when 397 and 866 nm fields are applied [l9[. 

The time dependence of the optical pulses is chosen as 
follows {h = 0, ii = (1 - a)T, t2=T and is = (2 - a)T 
where < a < 1): 



n'^it) = 




\^2{t)\' = 



■^'2. max 
■^'2, max 




{t < to) 
C0s'^TTt/2T) {to<t<t2) 
{t2<t< h) 

{t > ta), 

{t < to) 
{to<t< ti) 
cos"^ IT {t-ti)/2T (ti < i < ig) 

{t>t3), 



(3) 



where fi?,,^ and ili , 



(4) 



"max ""'" " '2 max ^^^ the maximum values for Q (t) 
and |5l2(t)p, respectively. In the experimental results 
given in this article, the values of a were empirically cho- 
sen to in between 0.55 and 0.85 so that effective popula- 
tion transfer was attained in experiments and/or numer- 
ical simulations. It is noted that the values of a used in 
this article are not thoroughly optimized ones but those 
that give relatively good performances among tested. 



IV. RESULTS FOR X- AND ^-ROTATIONS 

We show that arbitrary one-qubit operations can be re- 
alized with the geometric method introduced here. This 
can be accomplished by showing X- and Z-rotations in- 
dependently, because arbitrary operations can be decom- 
posed into sequences of such rotations. 

The procedure for A-rotations is as follows. The ion 
is first prepared in \u) by sideband cooling and optical 
pumping, then initialized to |0), and a STIRAP sequence 
is applied to perform a geometric gate. After this, the 
populations of the final states are analyzed by mapping 
the qubit states to the optical transition using a rect- 
angular TT pulse on |i) <->• \u) [i = 0,1,2). The popu- 
lations in 1 2) and \u) are also measured to verify that 



the gate operation is closed to the manifold spanned by 
the qubit states {|0) , |1)}. Oq and Oi are adjusted to be 
equivalent, and the peak values for (Oo,f^i,^2)/27r are 
~ (125, 125, 170) kHz. The total time for the STIRAP 
sequence is set to be ^120 /is by taking appropriate ac- 
count of the adiabaticity. a — 0.85 is used for the results 
given in this section. 

Figure ^a) shows the results for A-rotations. A si- 
nusoidal oscillation of population between |0) and |1) is 
observed. The sum of populations in the other states |2) 
and \u) is below 0.05 over the entire region, with which we 
confirm that the rotation operations are almost limited to 
the computational manifold spanned by {|0) , |1)}. The 
contrasts of the populations of |0) and |1) are obtained 
by fitting to 0.951 ±0.009 and 0.975 ±0.008, respectively. 

The procedure for observation of Z-rotations is similar 
to Ramsey interferometry. A superposition of the qubit 
states |0) and |1) is first prepared by application of a 7r/2 
pulse on |w)-|l) and a tt pulse on |m)-|0). A Z-rotation is 
then performed with STIRAP, which is followed by a tt 
pulse on |w)-|0) and a 7r/2 pulse on |u)-|l). The last pulse 
causes an interference between \u) and |1), which is de- 
tected by a projection measurement using fluorescence. 
Here the peak values for (Jig, Oi, il2)/27r are set to be 
-(200, 0, 170) kHz and the total time for the STIRAP se- 
quence is —120 /iS. Figure [^Jb) represents the results for 
Z-rotations, which shows a population oscillation against 
$ with an almost unit contrast. From a least-squares fit, 
the contrast is obtained to be 0.937 ± 0.010. 



V. EVALUATION OF FIDELITY 

In order to investigate the action of gate opera- 
tions more quantitatively, we also performed estima- 
tion of gate fidelities by using quantum state tomogra- 
phy. We performed fidelity estimation in the following 
three cases: (a) X-rotation by angle tt, (b) Z-rotation 
by angle tt, and (c) Hadamard gate [rotation around 
n — (l/\/2, 0, — 1/a/2) in the Bloch sphere by angle tt]. 
As explained before, general operations of the STIRAP 
gate can be described as the following unitary operator: 



C^STIRAp(^l 



,4>i,<^) 

cos — / + in(6'i, 0i) • (Tsin — 



(5) 



where n(0i,(/)i)= (sin 26*1 cos(/)i,— sin 26*1 sin 0i, cos 26*1). 
Using this general expression, the unitary operators for 
the three cases given above are respectively written as 
follows: 



Ux = t/sTiRAp(7r/4, 0,7r), 
Uz = C/sTiRAp(0,0,7r), 
Uh = C/sTiRAp(37r/8,0,7r). 



(6) 

(7) 
(8) 



The initial states in the three cases are respectively 
prepared as follows: (a)|i/'ox) = |0), {h)\tpoz) = (— « |0) + 



|l))/\/2, and (c)|7/'o_f/) = |0). In order to take into ac- 
count imperfect preparation of these initial states, we 
performed measurement of these states using the tech- 
nique of quantum state tomography. In the encoding 
adopted here, both the two states |0) and |1) in the com- 
putational subspace are in the same electronic state D^n- 
In order to discriminate these two states, we applied a 
mapping pulse between |0) and \u) before performing flu- 
orescence detection. The population out of the compu- 
tational manifold (population in |1) and |u)) before ap- 
plication of the mapping pulse, which is typically below 
0.05 as described before, is ignored here for simplicity. 
The mapping pulse itself produces a geometric phase rel- 
ative to the other states (|1) and |2)), which should be 
taken into account properly in the fidelity analysis. The 
mapping operation can be described as follows: 



RswA 



P,Om 



11 



|2)(2|. (9) 



This corresponds to a rotation by angle n around the x 
axis in the Bloch sphere of the two-level system {|0) , \u)}. 
The density matrix after application of the mapping 
pulse were reconstructed (using linear reconstruction) by 
fluorescence detection and an optional 7r/2 pulse prior 
to that. The populations along three orthogonal axes 
in the Bloch sphere, Pxii,Pyii and Pzu, are measured, 
where P^u {Pyii) is the population in -D5/2 after the 
mapping pulse and a 7r/2 pulse with the phase — 37r/2 
(0) on |u)~|l), and Pzu is the population in D5/2 state 
immediately after the mapping pulse. Using the Bloch 
vector r= {rx,ry,rz)= {l-2Pxu,l-2Pyu,l-2Pzii) the 
density operator after the mapping pulse is expressed as 

Pilots 

^ 1/ + ^(Iz.) (1| + |1) {u\) + 'f{-^ \u) (1| + ^ |1) {u\) 



+ ^{\u){u\-\l){l\), 



(10) 



and by using this the initial density operator is described 
as follows: 



- £.t 



Pi = RsWAP.OuPi.obsRsWAPflu- 



(11) 



In a similar way, the final density operator just after STI- 
RAP gate operations, py, is obtained in terms of popu- 
lations Pxif,Pyif and Pzif (defined in the same way as 
above) after a mapping pulse that follows the gate oper- 
ation. 

The gate fidelities in the three cases can be described 
using the density operators and the unitary operators, as 
follows: 



(12) 
(13) 

(14) 



Fx = tri^pfUxP^U}^ 
Fz = tr (^pfUzPtUlj , 

Fh = tT(^PfUHP^Ul 



The fidelities are explicitly written in terms of the initial 
and final populations, as follows: 

Fx = Pxli + Pxlf — 2PxliPxlf 

^Pyli ^ Pylf + '^PyliPylf 
+PzU+Pzlf-2PzuPzlf, (15) 

Fz — Pxli + Pxlf — 2PxliPxlf 

+Pyli + Pylf — 2PyiiPyif 
-PzU-Pzlf+2PzuPzlf, (16) 

Fh = Pxli + Pxlf — 2PxiiPxif 

+Pyli + Pzlf + 2PyiiPzlf 

+PzU + Pylf - 2PzuPylf - 1. (17) 

In order to determine the confidence intervals of the fi- 
delities, we obtained variances of the fidelities considering 
propagation of uncertainty based on the above expres- 
sion, as follows: 



V{Fx) 

= V{Pxu) + V{Pxlf) 

V{Pyu) + V{Py,j) - 

V{Pzu)+V{Pzlf)~ 

V{Fz) = V{Fx), 
V{Fh) 

= V{Pxu) + V{Pxlf)- 

V{Pyu) + V{Pzif) - 
V{Pzu) + V{Py,f) - 



^AP^,,V{Pxif)+AP^,fV{Pxu) 

AP^uV{Py,f) + AP^,fV{PyU) 

AP^,,V {Pzlf )+AP^,fV (Pzu), 



-AP^,,ViPxif)+4P^,fViPxu: 

■^PyuV{Pzlf)+^P^zlfV{Pyu) 

-AP^zuV{Pyif) + AP^^fV{Pzu) 



where V{...) represents the variance of each quantity. The 
variances of the populations [such as V{Pxii)] are simply 
determined here as the variances in binomial distribu- 
tions; e.g. V{Pxii) = Pxu{'^ - Pxii)lN, where N is the 
number of experiments per one measurement condition. 
The confidence intervals (68%) for the fidelities are deter- 
mined as the square roots of the variances of the fidelities. 

Table H] shows the results for the initial and final popu- 
lation measurements. The measurements were performed 
in essentially the same conditions as used in the previous 
section. Using the measured populations, the fidelities 
are estimated and shown in Table [TTl The second column 
in Table |TT] shows the fidelities reflecting both initial and 
flnal populations. We should note that these values are 
largely affected by imperfect initialization and analysis 
and therefore should be rather taken as the lower limits 
for the fldelities. 

For the purpose of reference, we also estimated the fl- 
delities when ideal initial or final states are assumed (the 
third and fourth columns in Table [TTl respectively). The 
fidelities for ideal initial states (the third column) repre- 
sent the goodness of the generated states without consid- 
ering initialization errors. These values that are better 
than those in the second column can be considered as 
the upper limits for the fidelities (when all operations 
other than the gates are assumed to be perfect). The 



Gate type P^u Pyu P^u Pxif 


Pylf Pzlf 


X-TT 0.510 0.549 0.026 0.463 


0.555 0.985 


Z-TT 0.509 0.045 0.502 0.472 


0.973 0.499 


Hadamard 0.497 0.515 0.025 0.485 


0.989 0.488 



TABLE I: Results of population measurements for fi- 
delity estimation. The initial {{Pxii,Pyii,Pzii}) and final 
{{Pxif , Pylf , Pzlf}) populations in the different three bases 
are measured for three different types of gate operations (see 
the text for details). Each result is obtained as the average of 
1500 experiments. 



Gate type 


Fidelity 


Fidelity Fidelity 
(ideal initial (ideal final 
states) states) 


X-TT 


0.965±0.038 


0.985±0.026 0.974±0.027 


Z-TT 


0.931±0.038 


0.973±0.036 0.955±0.039 


Hadamard 


0.965±0.038 


0.989±0.026 0.975±0.027 



TABLE II: Estimated fidelities. The confidence intervals are 
for 68% confidence level. The second column shows fidelities 
calculated from the measured initial and final populations. 
The third and fourth columns show estimated fidelities when 
ideal initial or final states are assumed, respectively, for the 
purpose of reference (see the text for details). 



fidelities for ideal final states (the fourth column) rep- 
resent the goodness of the initialization process. These 
values support the observation that the values in the sec- 
ond column is largely affected by imperfect initialization 
(and possibly imperfect analysis). 

We also performed fidelity analysis for STIRAP gates 
with rotation angle of 7r/2, which are more sensitive to 
coherence between STIRAP and square-pulse operations 
and hence to AC Stark shifts. Table ITTll shows the results 
for the initial and final population measurements for X 
and Z gates. Either | -ha;) = (|0)-|-|1))/V2 or |-|-z) = |0) as 
specified are used as the initial state. It should be noted 
that this time AC Stark shifts are not compensated at 
all, and relatively small values of the Rabi frequencies 
are used to suppress their effect. The peak values for 
(fio,^i,^2)/27r are set to be -(33.5,33.7,47.0) kHz for 
X gates and --(47.3,0,47.2) kHz for Z gates. The total 
time for the STIRAP sequence is —290 ^s and a = 0.75 is 
used. The results of population measurements are shown 
in Table IIIII Using these populations, the fidelities and 
their confidence intervals are estimated in the same way 
as above and shown in Table IIVI 

It should be noted that the fidelities in the case of ideal 
final states in Table HVl are higher than what arc given in 
Table HT] When taking the results in Table IIII[ we used 
a noise eater with a sample and hold capability to re- 
duce fluctuations in the amplitudes of the square pulses. 
We speculate that this helped improving the operational 
fidelity of the initialization. 



Gate type 
X-7r/2(|-|-a;) prep.; 
X-Tr/2{\+z) prep.; 
Z--k/2{\+x) prep.) 
Z-n/2{\+z) prep.) 



Pxli Pyli Pzli Pxlf Pylf Pzlf 

0.009 0.505 0.515 0.029 0.611 0.622 

0.526 0.540 0.006 0.595 0.916 0.363 

0.009 0.505 0.515 0.592 0.055 0.515 

0.526 0.540 0.006 0.546 0.534 0.041 



TABLE III: Results of population measurements for fidelity 
estimation of X and Z gates with it/2 rotations. Either 
\+x) = (jo) + ll))/\/2 or \+z) = |0) as specified are used 
as the initial state. Each result is obtained as the average of 
1500 experiments. Note that the same results are used as the 
initial populations for X and Z gates. 



Gate type Fidelity Fidelity Fidelity 

(ideal initial (ideal final 
states) states) 

X-7r/2(|-Ha;) prep.) 0.960±0.039 0.971±0.025 0.991±0.026 
X-tt/2{\+z) prep.) 0.905±0.037 0.916±0.026 0.994±0.026 
Z-7r/2(|+2:) prep.) 0.936±0.038 0.945±0.026 0.991±0.026 
Z-Tv/2{\+z) prep.) 0.952±0.038 0.959±0.026 0.994±0.026 

TABLE IV: Fidelities for gates with 7r/2 rotations estimated 
from measured populations in the previous table. 



VI. RESULTS FOR QUBIT ENCODING TO 
OPTICALLY SEPARATED LEVELS 



We have also attempted encoding of the qubit into op- 
tically separated levels. Use of optically separated lev- 
els may be advantageous when considering a combina- 
tion with optical two-qubit gate schemes @, [ifl, [ll| and 
use of less magnetic-field-sensitive transitions, such as 
Si/2{iTT'J = l/2)-i^5/2(™J = 1/2). This encoding uses 
an RF transition between the ground Zeeman sublevels 
[18j . Figure [Ije) shows the levels used for this encod- 
ing. {|0) , |1) , |2) , |u)} are encoded into Si/2{mj = 1/2), 
D^/ii'mj = -3/2), 05/2(111 J = -5/2) and 5i/2 (mj = 

— 1/2), respectively. 

Fig. [SJc) shows the results for X-rotations with qubit 
encoding to optically separated levels. Blue hollow cir- 
cles (red filled circles) represent the populations in |0) 
when |0) (|1)) is initially prepared. The contrasts are 
0.967 ± 0.012 and 0.916 ± 0.014 for the preparation of 
|0) and |1), respectively. Fig. [21[d) shows the result for 
Z-rotations, with a contrast of 0.886 ± 0.022. In these 
results, the peak values of {i}o,ili,fl2)/2TT are set to be 

- (130, 100, 100) kHz and - (130, 0, 100) kHz for X- and 
Z-rotations, respectively. The total time for the STIRAP 
sequence is approximately 128 fis. 

In ideal cases of X-rotations, the peak values of f^o and 
ill should be equal, which is not the case in the experi- 
ment described above. For a technical reason concerning 
the difference in the amplitude modulation of RF and 
optical fields, the temporal shapes of fto and fli are not 
proportional to each other in the present setup. The peak 
values given above are determined empirically so that 



they give correct X-rotations. This imperfection may be 
avoided by simply cahbrating the amplitude modulation 
process for either the RF or optical field so that JIq and 
rii give exactly the same time dependency. 



VII. DEMONSTRATION OF ROBUSTNESS 

The feature of the geometric phase gate is that it is 
expected to be robust against variations in the pulse area, 
namely the Rabi frequency and the pulse length. To 
study this feature, two experiments were conducted. One 
is the measurement of Jf-rotations with variation of the 
peak Rabi frequency and the illumination period. The 
other is the measurement of Z-rotations with variation 
of the peak Rabi frequency, while the illumination period 
is held constant. The ratios of the peak Rabi frequencies 
are held constant in each case. Both experiments are 
performed by measuring population oscillations as the 
rotation angle (/)2 is varied. 

The results for the former experiment is shown in Fig. 
[3] (a) and (b) for three values of the illumination period, 
with (a) fringe contrasts and (b) shifts. The ratios of the 
Rabi frequencies are held constant, and the peak values of 
ilo/27r are shown in the horizontal axes as representatives 
of the Rabi frequencies, a = 0.55 was used for the results 
given in Fig. [31 which gave relatively high contrasts for 
the case of 172/27r ^ 20 kHz. There is a region (above 
10^30 kHz) over which the absolute values of the fringe 
contrasts are almost independent of the Rabi frequency, 
which demonstrates the robustness. The fringe shifts are 
expected to start from zero and increase quadratically 
as the Rabi frequency is increased due to the increase of 
AC Stark shifts. The observed shifts shown in Fig. |3] (b) 
basically follow this expectation. 

The results for the latter experiment are shown in Fig. 
[3]Jc) and (d) with (c) fringe contrasts and (d) shifts. The 
fringe contrasts do not change appreciably as the peak 
Rabi frequency is varied in the region above 20 kHz. 

We have performed numerical simulations trying to ex- 
plain the loss of fringe contrasts and the shifts. Sim- 
ulations are performed using a Liouville equation with 
decay terms representing laser-frequency and magnetic- 
field fiucutuations. AC Stark shifts due to different 
Zeeman components and electric-dipole-allowed transi- 
tions {S1/2-P1/2, '§'1/2-/3/2 and -05/2-^3/2) are fully 
taken into account as time-dependent variation of the 
detunings. The results are shown as curves in Fig. |31 
where overall qualitative and partial quantitative agree- 
ments between experimental and numerical results are 
obtained. 



VIII. DISCUSSIONS 

The possible factors for loss of fidelity in the previ- 
ously described results includes laser frequency and mag- 
netic field fiuctuations. The laser linewidth in the current 



setup is ~ 300 Hz and the magnetic field fluctuations cor- 
responds to fluctuations of resonance frequencies of up to 
^ 30 Hz. The effects of thermal distribution of motional 
quantum numbers and intensity fluctuations are negligi- 
ble, because the scheme used here is robust against vari- 
ations in Rabi frequencies. 

We performed numerical simulations (similar to what 
is described in the previous section) to evaluate the fi- 
delity of some of the gates given in in f|Vl For X-rotation 
by angle n and Z-rotation by angle tt, the dominant infi- 
delity factor was the effect of magnetic field fluctuations 
(~1 %). Infidelity due to diabaticity was estimated to 
be ^--^0.2%. The contribution from laser-frequency fluc- 
tuations in this case was below 0.1 %. This can be ex- 
plained by the fact that the scheme is not sensitive to 
one-photon detunings but only to two-photon detunings 
(as will be described later in this section). The contri- 
bution from laser-frequency fluctuations may be larger 
when the encoding as in Wil is used. AC Stark shifts were 
compensated in this case by detuning one of the beams 
as described before. The possible fidelity loss from not 
performing compensation would amount to 20-30 % in 
our simulation. 

For the results of X- and Z-rotations by angle 7r/2, the 
largest infidelity factor was diabaticity (2-5 %), and the 
effect of magnetic field fluctuations was 1-3 %. Although 
no AC-Stark-shifts compensation was performed in this 
case, infidelity due to AC Stark shifts was as small as 
0-1 %, and the effect of laser-frequency fiuctuations was 
0.2-0.3 %. 

The scheme demonstrated here is useful not only for 
atomic systems (neutral or ionic), but also for other sys- 
tems including solid-state systems. The scheme is effec- 
tive for those systems where the coherence time is rel- 
atively long, but inhomogeneities of the excitation field 
intensities cannot be avoided, or for those where slow in- 
tensity fluctuations of excitation fields occur. It can also 
be used effectively in assuring equal operations in a large 
ensemble of particles. The requirements for the scheme 
are that a sufficient number of levels are available for im- 
plementation, and that energy shifts due to excitation 
fields (such as AC Stark shifts) can be avoided. 

Elimination of AC Stark shifts is not a straightforward 
task for the scheme presented here since a number of 
fields with time-dependent amplitudes are used. It can 
be confirmed from numerical simulations that the opera- 
tions of the STIRAP gates are not sensitive to one-photon 
detunings (i.e. detunings for |i)-|M) where i — 0,1 and 
2) but to every two-photon detuning, and hence all rel- 
ative shifts between two of |0), |1) and |2) should be 
minimized. This might be relatively complicated for e.g. 
^^Ca"*" since, in addition to adjacent Zeeman components 
in Si/2-D^/2, dipole-allowed transitions such as 6*1/2- 
Pi/2, S1/2-P3/2 and -D5/2-P3/2 also give rise to apprecia- 
ble AC Stark shifts t20|]. 

The scheme demonstrated here realizes single-qubit 
operations that are in general noncommutable to each 
other by using adiabatic manipulation of dark states. 



These amount to implementations of noncommutable 
(non-Abelian) holonomies, which can be a building block 
of HQC by Zanardi et al.^ and leads to application in 
other fields including neutral-atom systems |21j . 
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FIG. 1: (Color online) (a) Level scheme of a four- level sys- 
tem used to implement holonomic single-qubit operations, 
(b) Pulse sequence for holonomic single-qubit operations, (c) 
Level scheme used for ''"Ca"'" in the experiment, (d) Exper- 
imentally used pulse sequence, (e) Level scheme that imple- 
ments another four-level system containing one RF and two 
optical transitions, with qubit levels jO) , jl) separated by an 
optical frequency. 
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FIG. 2: (Color online) Results of single-qubit rotations (X- 
and Z-rotations). Variations of population are plotted against 
the phase shift in the middle of the gate pulse. The number 
of experiments per data points is 200 for experiments using 
three optical transitions [(a) and (b)] and 100 for those using 
one RF and two optical transitions [(c) and (d)]. The error 
bars represent errors in projection measurements, which are 
derived as standard deviations in binomial distributions with 
numbers of samples as just given, (a) Results for X-rotations 
using three optical transitions. Blue hollow circles, red filled 
circles, magenta crosses and black asterisks represent the pop- 
ulations in |0), |1), 1 2) and \u), respectively, (b) Results for 
Z-rotations using three optical transitions. The population 
in |0) is plotted, (c) Results for X-rotations using one RF 
and two optical transitions. Red filled circles (blue hollow 
circles) represent the populations in |1) and |2) when |0) (|1)) 
is initially prepared, (d) Results for Z-rotations using one RF 
and two optical transitions. The population in |1) and |2) is 
plotted. 
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FIG. 3: (Color online) Demonstration of the robustness of 
single-qubit operations. Each point represent the fringe con- 
trast and shift of the population oscillation obtained by vary- 
ing the rotation angle of a certain gate operation. Here the ro- 
tation angle is varied from to 2n in 17 steps, and the number 
of experiments per angle step is 50. The errors in the popula- 
tion measurements are estimated in the same way as in Fig. 
[2j and the fringe contrasts and shifts are obtained through 
weighted fits considering those errors. The error bars in the 
figure represents errors in the fitted parameters, a = 0.55 is 
used here. (a)Fringe contrasts and (b) shifts for X-rotations 
as the peak Rabi frequencies are varied. Here the ratios of 
the peak values are fixed. The horizontal axis represents the 
peak value of ^2 , and the blue circles, red triangles and green 
squares represent the results for pulse durations of 145, 290 
and 435 fis, respectively. Numerically calculated results for 
pulse durations of 145, 290 and 435 fis are also plotted as blue 
solid curves, red dashed curves and green dotted curves, re- 
spectively, (c) Fringe contrasts and (d) shifts for ^-rotations 
as the peak Rabi frequencies are varied. The ratios of the 
peak values are fixed. Numerically calculated results are also 
plotted as a blue solid curve. 



